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Abstract 

We show that the C 2 discrepancy of the explicitly constructed infi- 
nite sequences of points (xq, x\,x 2 , . . .) in [0, l) s over F2 introduced in 
[J. Dick, Walsh spaces containing smooth functions and quasi-Monte 
Carlo rules of arbitrary high order. SIAM J. Numer. Anal., 46, 1519- 
1553, 2008] satisfy 

C 2 , N ({x ,x 1 , . . < CsN- 1 (log N) s ' 2 for all N > 2, 

and 

C 2 ,2^({x ,x 1 ,. . . ,x 2 m^}) < C s 2" m m ( - S - 1 ^ 2 for all m > 1, 

where C s > is a constant independent of N and m. These results 
are best possible by lower bounds in [P.D. Proinov, On the L 2 dis- 
crepancy of some infinite sequences. Serdica, 11, 3-12, 1985] and [K. 
F. Roth, On irregularities of distribution. Mathematika, 1, 73-79, 
1954]. Further, for every N > 2 we explicitly construct finite point 
sets {y , . . . , f/jv-i) m [0> l) s such that 

£ 2 ,N({yo, Vi, ■ ■ ■ , Vn-x}) < CsN-HlogN)^ 2 . 

Another solution for finite point sets by a different construction was 
previously shown in [W. W. L. Chen and M. M. Skriganov, Explicit 
constructions in the classical mean squares problem in irregularity of 
point distribution. J. Reine Angew. Math., 545, 67-95, 2002]. 
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1 Introduction and Statement of the Main 
Results 

We study equidistribution properties of point sets in the s-dimensional unit- 
cube [0, l) s measured by their C 2 discrepancy (see [21 [121 EES [2Ql 121]). For 
a finite set Vn, s = { x o, • • • , x n-i} of points in the s- dimensional unit-cube 
[0, l) s the local discrepancy function is defined as 

a / x A N ([0,u>),V N , s ) 

A(wi, ...,uj s ) = — wi • • -u s , 

where u> = (ui, . . . ,u s ) 6 [0, l] s and An([0,lj),Vn,s) denotes the number 
of indices n with x n 6 [0, u-\) x ■■• x [0,u; s ) =: [0, a>). The discrepancy 
function measures the difference of the portion of points in an axis parallel 
box containing the origin and the volume of this box. Hence it is a measure 
of the irregularity of distribution of a point set in [0, l) s . 
The C 2 discrepancy of Vn, s is defined as 

Un(Pn,s)= (I |A(u,)| 2 d^ ' . (1) 



'[0,1]- 

For an infinite sequence S s in [0, l) s the C 2 discrepancy C2,n(<S s ) is the C 2 
discrepancy of the first iV elements of S s . It is well known that a sequence 
is uniformly distributed modulo one if and only if its C 2 discrepancy tends 
to zero for growing N. Furthermore, the C 2 discrepancy can also be linked 
to the integration error of a quasi-Monte Carlo rule, see, e.g. [TJl 1251 ES] for 
the error in the worst-case setting and [39J for the average case setting. 

A lower bound on the C 2 discrepancy of finite point sets has been shown 
by Roth [32] which states that for any s G N (the set of positive integers) 
there exists a number c s > depending only on s, such that for every point 
set Vn,s in [0, l) s consisting of iV > 2 points we have 

r(T> ^ (log AO^- 1 )/ 2 

C 2 {V N , S ) > c s — . (2) 
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This lower bound is best possible in the order of magnitude in N as shown 
first by Davenport [7] for s = 2 and then by Roth [331 EI] and Frolov [T7] 
for arbitrary dimensions sGN. Davenport used point sets consisting of the 
2N elements ({±na}, n/N) for 1 < n < N, where JVeN and a has a con- 
tinued fraction expansion with bounded partial quotients. Further examples 
of two-dimensional point sets with best possible order of £ 2 discrepancy can 

be found in [H E2 CESJ [2U ED]- 0n the other hand ' Roth ' s Proof 
for dimensions s > 2 is a pure existence result obtained by averaging argu- 
ments. Explicit constructions of point sets achieving the best possible order 
of convergence have been a longstanding open problem. Finally, a solution 
was given by Chen and Skriganov [I] who, for every integer N > 2 and every 
dimension s G N, gave for the first time explicit constructions of finite point 
sets consisting of N points in [0, l) s whose £2 discrepancy achieves an order 
of convergence of (log N)^* 1 " 2 /N. Their construction uses a finite field ¥ p 
of order p with p > 2s 2 . We also refer to [5] where the arguments from [I] are 
considerably simplified and to the overview in [121 Chapter 16]. The result 
in [1] was extended to the C p discrepancy for 1 < p < 00 by Skriganov [35J. 

On the other hand, it was shown by Proinov [21] that for an infinite 
sequence S s = (xq, Xi, . . .) of points in [0, l) s there is a constant c' s > such 
that 

£2,MW asi, . . . , a?jv_i}) > c s — 

for infinitely many values of N. This lower bound is known to be best possible 
in dimension s — 1. One-dimensional infinite sequences whose £2 discrep- 
ancy satisfy a bound of order y/log N /N for every N > 2 were given in, e.g. 
[31 ESI Ell ESI ED]- These constructions are mainly based on the symmetriza- 
tion of sequences (also called reflection principle). It is widely believed that 
Proinov's lower bound is also best possible for arbitrary dimensions s, which 
we show in this paper. The authors are not aware of any other construction 
which matches Proinov's lower bound. 



1.1 The Main Results 

In this paper we prove two main results: We provide explicit constructions 
of infinite sequences in [0, l) s for which the first N > 2 points achieve a £2 
discrepancy of order (log N) s ^ 2 /N for arbitrary s G N. This result is best 
possible by the lower bound of Proinov [29] . 
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Furthermore, for any integer N > 2 and any dimension s G N, we give an 
explicit construction of a finite point set of N elements in the s-dimensional 
unit cube with the optimal rate of convergence for the L 2 discrepancy in the 
sense of the lower bound of Roth. Our construction is completely different 
from the construction of Chen and Skriganov [I]. In contrast to [I] where the 
construction uses a finite field ¥ p with p > 2s 2 our method is, independent 
of the dimension s, based on the finite field F 2 of order two. Furthermore, 
our result does not use the Davenport reflection principle [7] and also does 
not use the 'self-averaging' property from [4]. Instead it is based on higher 
order digital nets and sequences from [8j [9]. 

In our proofs we do not keep track of constants which depend only on the 
dimension s since they are significantly larger than the constants obtained 
in [11]. Therefore, in the following, we write A(N,s) <C S B(N,s) if there 
is a constant c s > which depends only on s (and not on iV or m through 
N = 2 m ) such that A(N, s) < c s B(N, s). 

Theorem 1 For any sGN one can explicitly construct an infinite sequence 
S s = (x , x±, x 2 , . . .) of points in [0, l) s such that for all N > 2 we have 

(log N)^ 1 ^ 2 , (loeiVW 2 

C 2 , N ({xo, a*-i}) « s VS(Nj « s { g N } , 

where S(N) is the sum- of- digits function of N in base 2 representation, i.e. 
if N = 2 ni +2 n2 + ■■■ + 2 nr with m > n 2 > ■ ■ ■ > n r > 0, then S(N) = r. 
Obviously, we have S{N) < 1 + (log N)/ (log 2) for all N G N. 

Remark 1 It follows from (22J, Corollary 3] that for any e > we have 



lim — 

M-s-oo M 



log M log M 

0<N<M : <S{N) < + * 



1. 



21og2 v ' v '2^2 

Hence the density of N G N for which S(N) is at least of order logiV is equal 
to one. More precise results on the distribution of the sum-of-digits function 
can be obtained, e.g., from [H 123] . 

The above construction can also be used to obtain the following result 
for finite point sets, which was first shown in [1] by a different construction. 

Corollary 1 For any s G N and any integer N > 2 one can explicitly con- 
struct a point set Vn,s consisting of N elements in [0, l) s such that 

(log AQ(-i)/ 2 

i->2,N(,rN,s) • 
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1.2 Explicit constructions of sequences and point sets 

We now present explicit constructions of sequences and point sets satisfying 
Theorem [1] and Corollary [TJ 

The construction of sequences S s = (xq, X\, . . .) in [0, l) s satisfying The- 
orem [TJ was introduced in [HIE] and is based on linear algebra over the finite 
field F 2 of order 2 (we identify F 2 with the set {0, 1} equipped with the 
arithmetic operations modulo 2). 

Let Ca,...,C s G F 2 xN beNxN matrices over F 2 . Let n = n + 
n\2 + • • • + n m _i2 m_1 G No be the binary expansion of n and set n = 
(n , rii, . . . , n m _i, 0, 0, . . .) T G F 2 . Then define 

Xj >n = Cjfi for j — 1, ... , s. 

Let Xj >n = (xj )ri: i, Xj : n )2 , ■ ■ -) T an d define 

Then the nth point x n of the sequence S s is given by x n = (xi jn , . . . , x S)n ). 

Explicit constructions of suitable generating matrices Ci,...,C s where 
obtained by Sobol' |37J . Niederreiter [26], Niederreiter-Xing [28J and others 
(see also [T21 Chapter 8]). We briefly describe a special case of Sobol's and 
Niederreiter's construction of the generating matrices. Let Cj = (cj^/)k/>i 
with Cjfci G F 2 . Let pi = x and pi G ¥ 2 [x] be the {i — l)th primitive poly- 
nomial in a list of primitive polynomials over F 2 that is sorted in increasing 
order according to their degree ej = deg(pj), that is, t\ < e 2 < • • • < e s (the 
ordering of polynomials with the same degree is irrelevant). Consider the 
Laurent series expansion 

— — = ^a £ (z, j, z)x- e G F 2 ((a;- 1 )). 

Then for 1 < j < s, k, i > 1 we define 

Cj,k,e = ae{i,j, z), 

where k — 1 = (i — l)ej + z with integers v and i satisfying < z < ej. 

Using the generating matrices defined this way, we obtain an infinite 
sequence S s = (x , x%, . . .) for any dimension s G N. 

To obtain the sequence which satisfies Theorem [TJ we need the following 
definition. 
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Definition 1 For d G N the digit interlacing composition (with interlacing 
factor d) is defined by 

® d : [0,l) d -> [0,1) 

oo d 
a=l r=l 

where x r = £ rjl 2 _1 + ^ r 2 2~ 2 + ■ ■ • for 1 < r < d. We also define this function 
for vectors by setting 

$} d :[0,l) ds -> [0,l) s 

(xi, . . . ,Xds) ^ (^(xi,...,x d ),...,^(a;( s _i) d+ i,...,x sd )), 
for point sets Vnm = {xq, x%, . . . , x^-i} C [0, l) sd by 

0d(7V,«i) = {®d(Xo), • • • , ^(*iV-l)} C [0, l) s 

and sequences iS s d = (x , x±, . . .) with a; n G [0, by 

Then the sequence ^(Ssa) satisfies the bounds in Theorem [TJ More 
generally, any sequence 3td(<Sds) with d > 5 satisfies the bound. 

Note that the interlacing can also be applied to the generating matri- 
ces Ci, . . . , C s directly as described in P Section 4.4], thus one can obtain 
generating matrices E\,...,E S G F2 XN which generate a digital sequence 
satisfying Theorem [TJ 

To construct finite point sets for any integer N > 2 we proceed the 
following way. Let m G N be such that N < 2 m and let x Q , X\, . . . , x?m-\ G 
[0, l] 3s_1 be the first 2 m points from the Sobol or Niederreiter sequence in 
dimension 3s — 1 as introduced above with p\ — x and P2 = 1 + x. Let x n = 
(x ltn , . . . , x 3s _i, n ) and define y n = (n2~ m , x 1>n , . . . , £ 3 s-i,n) G [0, l) 3s . Let 
now V2 m ,s = {^3(1/0), ^(Vi), ■ ■ ■ 1 ^3(2/2™)}- To obtain a point set consisting 
of N points we use a propagation rule introduced in |4] (see also [HI p. 512]): 
The subset 



Vn,s '■— 7V\s n 



contains exactly N points. Then we define the point set 



rn 



P \ ..- := { ( jjrXi-->-2 '•, } : (-''[- -r 2 r.J G V \..-\ ■ (3) 
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Then Vn,s satisfies the bound in Corollary [TJ We remark that Chen and 
Skriganov [4 J applied the same propagation rule but to a different point set. 

1.3 The general construction principle 

The constructions of [Bl E] were introduced to obtain explicit examples of 
so-called higher order digital sequences. In fact, our proof shows that any 
higher order digital sequence [U Definition 4.3 and Definition 4.8] of order at 
least 5 satisfies Theorem [TJ We state simplified versions of these definitions 
here which are sufficient for our purpose. 

Definition 2 Let m, a G N and let t be an integer such that < t < am. 
Let d,...,C a e F™ xm with Cj = (c jA , . . . , c iim ) T , i.e., c jti G F™ is the zth 
row vector of the matrix Cj. If for all 1 < ij >v . < ■ ■ ■ < ij\ < m with 

s mm(uj,a) 

Yl hi <ocm-t 
j=i i=i 

the vectors 

are linearly independent over F 2 , then the digital net with generating matrices 
Ci, . . . , C s is called an order a digital (t, m, s)-net over F 2 . 

Definition 3 Let a G N and let t > be an integer. Let d, . . . , C s G F 2 xN 
and let Cj )TnXm denote the left upper mxm submatrix of Cj. If for all m > t 
the matrices Ci imxm , . . . , C s>mXm generate an order a digital (t, m, s)-net over 
F 2 , then the digital sequence with generating matrices Ci, . . . , C s is called an 
order a digital (t, s)- sequence over¥ 2 . 

It follows from the results in [HI Theorem 4.11 and Theorem 4.12] for the 
finite field F 2 that the sequence @d(Sds) is an order d digital (t, s)-sequence 
(which shows that t is independent of m). 

The proof of Theorem [TJ shows that every order a digital (t, s)-sequence 
over F 2 with a > 5 satisfies Theorem[TJ We conjecture that a > 2 is sufficient. 
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2 Walsh series representation of the squared 
£2 discrepancy 

2.1 Walsh functions 

We introduce Walsh functions in base 2 (see [6l EE]), which will be the main 
tool in our analysis of the £2 discrepancy. In the following let No denote the 
set of non-negative integers. For k G No with base 2 representation 

k = Ka^- 1 + --- + k 1 2 + k , 

with Kj G {0, 1}, we define the Walsh function wal^ : [0, 1) — > { — 1, 1} by 

wal fc (x) := (_i)*i"o+-+*«*.-i j 

for x G [0, 1) with base 2 representation x = + §§ + •• • (unique in the 
sense that infinitely many of the x^ must be zero). 

For dimension s > 2, xi, . . . , x s G [0, 1) and k±, . . . , k s G No we define 
wal fcl ,..., fcs : [0,l) s ^{-l,l}by 

s 

wal fcli ... ifcs (xi, . . . ,x s ) := JJwal fci (a: i ). 

For vectors k = (ki, . . . , k s ) G Nq and x = (xi, . . . , x s ) G [0, l) s we write 

wal fc (aj) := wal fcli ... ifcs (xi, . . .,x s ). 

Further properties of Walsh functions can be found in [12], Appendix A]. 

We introduce some notation. By © we denote the digit-wise addition 
modulo 2 , i.e., for x = % and V = w we P ut 

00 

x @ y '■= 7 — , where Zi := X{ + yi (mod 2). 

Further we call x G [0, 1) a dyadic rational if it can be written in a finite base 
2 expansion. For vectors x,y G [0, l) s we set a? © y = (xi © yi, . . . ,x s © y s ). 

Let {#0, • • • , a?2 m -i} in [0, l) s be a digital net over F 2 and let cr G [0, l) s . 
Then we call the point set {x © cr, . . . . X'2 m — 1 © <r} a digitally shifted digital 
net over F 2 . 
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2.2 The Walsh series expansion of the £2 discrepancy 

The squared £2 discrepancy of a point set Vn,s = {xo, ■ ■ ■ ,&n-i} can be 
viewed as a function of {xq, • • • , xn_i}, i.e. a function of Ns variables: 



In the following we give the Walsh series expansion of the function C\ N . 

There is a well known formula for the £2 discrepancy of a point set by 
Warnock (see, for example, [T2| Proposition 2.15]), which we use in the fol- 
lowing to obtain the Walsh series expansion of C-2,n(Pn,s)- This formula is 
given in the following proposition. 

Proposition 1 Let Vn,s = {x®, . . . , a?jv-i} be a point set in [0, l) s . Then we 
have 



where x n j is the j-th component of the point x n . 

The following result deals with the Walsh series expansion of the indicator 
function l[o |£i ,)(x) (which is 1 for < x < cu and otherwise) and is a conse- 
quence of [211 Lemma 2 and Lemma 3] (see also [T2l Proof of Lemma 3.29]). 



Here and in the following let k > have base 2 representation k = 
k + k x 2 H h /t a -22 a ~ 2 + 2 a ~ 1 with m G {0, 1}. We define the weight 



^■2,nCPn,s) — ^^({^Cb • • • > x N _i}). 




function 




For any < x < 1 we have 




wal 2 M(*)-i (x) 



— rrrwal 
2 r+2 



r=l 




) 
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As the above formula gives the k-th Walsh coefficient of l[ 0jW )(a;) as a 
function of u we obtain the Walsh series expansion 

Mo,oj)(x) ^ l-x 

°° 1 ( °° 1 \ 
+ 2 M(fc)+i S ^wal t+2r+( ,w-i(i) - wal fce2fl(fe )-i(x) wal fc (w). 

fc=l \r=l / 

The well-known Parseval identity yields the equality 

/ l [0 ,u,)(^)l[o, w )(y) dw 
Jo 

00 1 ( 00 1 \ 

= (1 - - y) + XI 2 2Mfc)+2 walfce2M(*)-i(a:) - Yl ^ wal fc+r+"W-iW 1 

fc=l V r=l / 

< 00 ^ N 

wal fce2M(fc )-i(y) - ^ ^wal fe+2 , +M( fe)-i(y) 



' 2 T 

r=l 



Now we consider the £ 2 discrepancy of a point set 7-V,s = {xq, ■ ■ ■ , Xn-i} 

in [0, l) s . Using the fact that A N ([0, u), V N , S ) = T,n=o IlJ=i Mow)( x "j) [t 
follows that 



A N ([0,u),V N ," 2 



Combining the last two equations we obtain 
fA N ([0,oj),V NyS y 2 

[0,l] a V 



j JV-i s „! 

dW = ^ E 11/ l[0 W )(^)l[0^)(^mj) dWj. 
n,m=0 i=l ^ 



\ - I s 

1 ~~ ;r n,j)(l — x m,j) 



n,m=0 j=l 

°° 1 / " \ 

' S 2 2M(fc)+2 ( wal fee2M(fc)-i(^nj) - X ^wal fc+2 r +MW -i(a; nii ) 

fc=l V r=l 

/ 00 1 

x I wal fce2M ( fc )-i(x m j) - ^ — wal fc+2 r+ M (fc)-i(x m j) 



1 27 

r=l 



(4) 
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Note that we have J [Q l]s ( Ajv([0 '^ ) ' 7 ^ s) ) = J2n,m=o ITU min(l-x njj , 1- 
x m,j), which follows from the proof of Warnock's formula for the £ 2 discrep- 
ancy. 

The Walsh series representation of (1 — x n j){l — x m j) can easily be found. 
For example it was shown in (T2J Lemma A. 22] that 



1 

x - - = - 7^j-wal 2 a-i (x) . (5) 

a=l 



Using (J4]) together with the last equality we obtain the Walsh series repre- 

sentation of / [0 1]s ( A "«°-ff-^.'> ) 2 du>. 

Using again ([5]) and Proposition [T] we can now obtain the Walsh series 
expansion of the squared £ 2 discrepancy, which is given by 



1 2 N ~ 1 s ( I °° 1 
3l-jv£Il 3+E ^ wal 2»" 1 

n=0 j=l \ a=l 



£ ^^ Wal 2«-i+2«'-i( X n,i) ) 
l<a<a' / 



1 N ~ l S \ (l °° 1 \ / 1 °° 1 

+ AT2- Ell ( 2 + S 2^T wal2a - l(x "^ ) ( 2 + S 2^+1^-1 (s m>i ) 

n,m=0 j'=l L \ o=l / \ a=l 

00 ^ / 00 ^ 

+ 2 2 M(fc)+2 ( wal fee2M(fe)-i(^n,i) - ^Walfc+jjr+MW-i^nj) 
fc=l \ r=l 

x ^wal fcffi2M( fc)-i(x m) j) — ^ ^r w al fc+2 r+ M (fe)-i(a;mj) 

The following lemma can now be obtained upon comparing coefficients. 
Lemma 1 For any Vn, s = { x o, • • • > x n-i} in [0, l) s we obtain 

N-l 

3 s N 



1 2 N ^ 

£InCPn, s ) = — - Yl Yl r ( fc ' °) wal fc(a;„) 



n=0 ke% 
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1 N ^ 

+ Jp KM)wal fc (# n )wah(cc r 



n,m=0 feJeNS 



where k = (ki, . . . , A; s ) ; / = (Zi, . . . , r(fc, Z) = n^=i r (%> (?')■ Further we 
have r(k, I) = r(l, k) and for non-negative integers < I < k with k = 



2" 1 - 1 + • • • + 2 a "~ 1 with at > 
b\ > ■ ■ ■ > b v > we have 



> a v > and I = 2 hl ~ l + • • • + 2 



b w — 1 



with 



( l 

3 



r(k,l) 



i 

2°i+ 2 



ifk = I = 0, 

if v = 1 and I = 0, 



2 ai +a 2 +2 ifv = 2 and I = 0, 
^+^2+2 ifv = w + 2>2anda 3 = b 1 , 



ifk = l>0, 
: i i- ifv = w,a 1 ^b 1 anda 2 = b 2 , 



3-4 a i 



otherwise. 

The above formula can be simplified further. 

Lemma 2 27ie squared L 2 discrepancy of a point set Vn,s 
in [0, l) s can fre written as 



> flu — 0^-2) 



JV-1 



iV-1 



AT.sJ 



fc,zeNg\{o} 



n=0 



m=0 



where the coefficients r(k,l) are given as in LemmaUl 
If T y 2 m , s is a digital net over ¥ 2 we have 

C 2 22m (V 2 m :S ) = ^ r (M), 

where T>* =T>\ {0} and where T> is the so-called dual net given by 

V = {(&!, . . . , s) e Ng : Cjh + ■■■ + Cjk s = 0}, 
where for fc G N„ with base 2 expansion k = k + k x 2 + k 2 2 2 + ■ ■ • we pit£ 

IfV 2 m s is a digital net over F 2 digitally shifted by a e [0, l) s we /mve 



£ 2 

*-'2,2' 7 



(7Vs«) = 2j ?"(M)wal fe (cr)waL(< 



where T>* denotes the dual net excluding 0. 
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Proof. We have 



* n=0 fceNg 

+ 7\T2 r (M)wal fc (a; n )waL.(;E m ) 

n,m=0 fe,JeNg 

^ N-l Af-1 
fc,i£Ng\{0} n=0 m=0 

from which the first part follows. The second part follows immediately from 
the character property of digital nets and the above equality. The third part 
follows in the same manner as the second part using the additional equality 
walfc(a; © a) = walfc(ai)walfc(cr). □ 

3 The proof of Theorem U 

We give the proof of our main result. 

Proof. Let S s = (x , x%, . . .) in [0, l) s be an order a digital sequence over 
F 2 with a > 5. By [HI Theorem 3.3 (iii)] it follows that the sequence is 
an order a' digital (£(«'), s)-sequence over F 2 for 1 < a' < a. Let t = 
max{t(l), t(2), t(3), t(4), t(5)}, that is, the sequence is an order a' digital 
(t, s) sequence over F 2 for all 1 < a' < 5. Note that the sequence explicitly 
constructed in Section 11.21 is an order 5 digital (t, s)-sequence over F 2 and 
therefore the proof applies to this case. Further note that the first 2 m points 
of this sequence are an order a' digital (t, m, s)-net over F 2 for 1 < a' < 5 
by Definition [3l 

We use the first part of Lemma [2] to obtain 

. N-l ^ N-l 

£1,n(Ss)= Yl r(k,l)—^w&\ k (x n )—^2wali(x m ). 

k,l£W Q \{0} n=0 m=0 

Let N = 2 mi +2 m2 + ■■■ + 2 mr with m 1 > m 2 > ■ ■ ■ > m r > (hence 
r = S(N)). Then the points 

a; 2 m 1 _ ) — |_ 2 m i-i , • • • , x_i + 2 m i-\ — \-2 m i j 
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where for i = 1 we define 2 mi + • • • + 2 m '~ 1 = 0, are a digitally shifted 
order a' digital (t, mj,s)-net over F 2 for 1 < a/ < 5 (see Definition [3] and 
[T2"j Section 4.4.6]). The generating matrices Ci >miXmi , . . . , C SjTOiXmj of this 
digital net are the left upper rrii x m 8 submatrices of the generating matrices 
Ci, . . . , C s of the digital sequence. Assume that the digital shift is given by 
cTj. Let T>i denote the dual net corresponding to this digital net, i.e., 

Di = {k = (&i, ...,/:,) eNJ : C lm . xm .ki + • • • + C srn . Xm .k s = 0}, 

where for k G No with base 2 expansion k = n + k%2 + k,2^ 2 + • • • we set 
= (ko, K\, . . . , /t mi _i) T . Set £>* = T>i\ {0}. Then we have 

1 JV-l r 1 -l+2-i +... +2 »N 

-^wal fe (a; n ) — ^ wal fc (aj n ), 

n=0 i=l n=2 m iH \-2 nli - 1 

where again for i — 1 we set 2 mi + ■ ■ • + 2 mi - 1 = 0, and by [121 Lemma 4.75] 

1 -l+2 m i+---+2 m i ( . . 

1 v - i / \ J wal fc (crj) it fe e P i; 

— ^ wal*(«„) = | Q if fe ^ 

n=2 m lH |-2 m »-i k 

Therefore 

^ AT-l - JV-1 

^2,7v(^) = 5^ r ( fc ' Z )^y 5^ ^fe^n)^ (^m) 
fc,ZeNg\{0} n=0 m=0 



— — ^ r(fc,Z)wal fc (cr i )wali(cr i 



N N 
i,i'=i kev*,iev*, 



< v — — — V irffc n 



iV 

i,i'=l keV'ltV* 

We introduce some notation. Let 

J i>v = {{k,l) eV* xV*, : r{k,l) ^U} 

and 

Ji t i>(z) = {(fe, Z) G ,7 M / : //(fe) + fi(l) = z}. 

Recall that the function \i : No — > No is defined by /i(0) = and for k 
k + k x 2 H h K a _ 2 2 a ~ 2 + 2"- 1 with G {0, 1} by n{k) = a. 
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For (k, I) G J it i> we have 

\r(k,l)\ < 

Thus we have 



3s2M(fe)+M(«) ' 



AT N ^ 2^( fe )+M0' 

i,i'=i (fc,i)eJi,i/ 

It follows from k G D* that > rrii — t and from Z G "£>*/ it follows also 
that > mj' — £ and hence ji{k) + /x(Z) > + r?ty — 2£ + 2. 
Thus we have 

E l = \Ji,i'( z )\ 
2M(fc)+/i(0 " 2 2 

(k,l)eJ iA t z=m i +m i ,-2i+2 

Let fcj, G No- In the following we simultaneously use two different 
notations for the binary expansion of kj and lj. First let 

kj = 2 < ^ 1_1 + --- + 2 a ^~ 1 

with aj ; i > ■ • • > a^. > and 

Z j = 2^ 1 - 1 + --- + 2 6 ^ -1 

with bj t i > ■ ■ ■ > bj^j > 0. Thus Vj denotes the number of nonzero digits of 
kj and uij denotes the number of nonzero digits of lj. For kj = we use the 
convention that Vj = and a^i = 0. Further we set a^.+j = bj^j+i = for 
i>0. 

Further we also use the notation 

/•> • hi 2 ^" ; 

with binary digits fc^j G {0, 1}. Thus 

^ f 1 if i — dj tV for some 1 < v < Vj, 
j ' 1 \ otherwise. 

Analogously we write 

h = hfl + ^',i 2 1 - O-'',.! i-'' 1 1 
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with binary digits G {0, 1}. Thus 

^ f 1 if i = for some 1 < u> < Wj, 
J ' 1 1 otherwise. 

Recall that r(k,l) = Yij=i r (^j^j)- For r(kj,lj) ^ it follows that in 
some sense kj and /.,• cannot be too different. Let us elaborate in more detail: 
Assume that r(kj,lj) ^ 0. Now Lemma [T] implies that in order for r(kj,lj) 
not to be we must have < \vj — Wj\ < 2. Further we must have: 

(z) \vj - Wj\ = a jj2 = b jj2 , • • • , a^ 3 = b j>v . , 
(ii) \vj — Wj \ = 1 =>- kj = or = 0, 

(Hi) \vj — Wj \ = 2 =>> if = w j + 2 then 0^3 = bj t i, . . . , o,-^ = bj^ j 

if = Vj + 2 then fy i3 = a jt 1, . . . , 6^ = a^.. 

If \vj — Wj\ > 2 we always have r(kj, lj) = 0. 

For given (k, I) G J%^(z) we define the following sets for —2 < r < 2: 

«r = {j G {1, . . . , s} : vj = Wj + r}. 

Note that a T (1 a T > =0 for r /t' and (Jr=-2 ^ = • • • > s ) by Lemma [TJ 
Then we have 

1. For j G «2 we have lj ti = kj ti for < i < dj i2 — 1; 

2. For j G ai we have J 3 - = and kj = 2 ajVi ~ 1 ; 

3. For j G ao we have kj ti = lj ti for < i < minjaj^, 6^1} — 1; 

4. For j G a_i we have lj = 2 6j > 1 ~ 1 and kj = 0; 

5. For j G a_2 we have Ij^ = kjj for < i < bj 2 — 1. 

Thus in all cases we have fc^j = ^ for < i < minla^ — 1, bj j2 — 1}. We set 
now 

hj t i = kjj = lj t i for all 1 < j < s and < % < minja-,-^ — 1, bj >2 — 1}, 
and for Uj = mm{a,j t2 — 1, bj t2 — 1} we set 

hj = h jfl + h jtl 2 + ■■■ + hj ]Uj „x2 u i~ l for 1 < j < s 
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if Uj > and hj = otherwise. Thus we only need to consider the cases 
where 

kj = hj + [2 a ^~ 1 \ + [2^ 1 - 1 J, 
\. = hj + [2 6 ^ 2 - 1 J + L2^ 1 ~ 1 J 

for 1 < j < s. 

We now study the number of elements in the set Ji^{z). Assume that 
i < i' . Assume that (k, I) G Ji,i'(z), z = z\ + z 2 with = Z\ and /i(Z) = z 2 . 
Let 

R(k,z 2 ) = {IE V* : (k,l) E J LV {z)^{l) = z 2 }. 

If < mi — t or z 2 < — t then z 2 ) = 0, thus we only need to 
consider the case where > rrii — t and ^2 > ^i' — t. Thus z = zi + z 2 > 
rrii + nti' — 2t + 2. We have 



z— m,-/+i— 1 

\\keV* : (x(fc) = zi}| 

fcex>* 

zi=nii—t+l ,, % 1 



Um'(*)I< £ : /i(fe) = zi}| max \R{k,z-z{)\. (6) 



We now study the two factors from the above sum separately in two steps: 

Step I The number of k = (ki, . . . , k s ) £ T>* with fi(k) = z\ has been 
studied in [TTj . Assume first that fcj > for 1 < j < s. The case where 
one or more of the kj's are zero follows by the same arguments. Let 
£(Vl, v 2 , . . . ,v s ) denote the number of such k = (hi, . . . , k s ) 6 V* with 
/x(fcj) = a>j,i = Vj. Then fc e P* implies that 

ci,i/ii, H h ci 

,ai,2-l^l,ai, 2 -2 + c l,ai,2 + Ci )fltl x + 
C 2; \h 2 fl + • ■ ■ + C2 >a2i2 _i/l2,a 2 ,2-2 + C2,a 2 ,2 + c 2,a 2 ,i + 

: ^ (7) 

Cs,ih s fl ~ ~ ~ ~\~ c S CLg 2—\h s as 2— 2 ~\~ c S CLg 2 -|- c SCLs ^ 0, 

where c^j G F™* denotes the ith row vector of the matrix Cj jmiXmi . 
Since by the digital (t,rrii, s)-net property the vectors 

■ ■ ■ j ■ ■ ■ j C S> 1, ■ ■ ■ j C S)? ; 3 

are linearly independent as long as v\ + ■ ■ ■ + v s < — t, we must have 
f i H h t> s > m, - t + 1. (8) 
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Let now A denote the rrii X ((^i — 1) + • • • + (v a — 1)) matrix with column 
vectors cj tl , . . . , cj^, . . . , cj 1? . . . , cj^, i.e., 

^ := ( c l,l; • • • j c 1 )Vi -1j • • • ) c s,l) • • • 5 C s,v a -l)- 



Further let 

f--=Sl Vl + --- + cl Vs e¥" 2 

and 

k '■= (klfl, • • • , fel,t; 1 -2) • • • j ^s,0; • • • 5 ^s,-u s -2) T £ ^ 



(«i-l)+-+(v 8 -l) 
2 



Then the linear system of equations (J7J) can be written as 

Ak = f (9) 

and hence 

E( Vl , ...,«.)= J] 1 = |{fc G F (^-i)+~+<«.-i) . A fc = / } |. 



S6F («i-D+-+(«.-i) 



By the definition of the matrix A and since Ci m , Xm ., . . . , C s m^xm^ 
the generating matrices of a digital (£, m;, s)-net over F2 we have 



rank^A) 



(vi - 1) + • • • + (v. - 1) if (vi - 1) + • • • + (v s - 1) < mi - £, 
> rrij — t otherwise. 



Let L denote the linear space of solutions of the homogeneous system 
Ak = and let dim(L) denote the dimension of L. Then it follows that 



dim(L) 



if V\ + • — h v s < rrii — t + s, 

< v% + • • • + v s — rrii + t — s otherwise. 



Hence if v\ + • • • + v s < rrii — t + s we find that the system ([9]) has at 
most 1 solution and if v 1 + • • • + v s > rrii — t + s the system (jHJ) has at 
most 2 Vl+ "' +Vs ~ mi+t ~ s solutions, i.e., 



E(vi, ...,v s ) < 



1 if Vi H — • + v s < rrii — t + s, 

2 v 1+ ...+v a -m t +t- S iivi + ... + Vs>mi _ t + s ^ 
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In the following let denote the binomial coefficient, where we set 
(™) = if k > n. Thus we have 

\{keV* : k J >0,l<j<s,fi(k) = z 1 }\ 

/ (^T 1 ) ifz 1 <m,-t + S , 

\ (^+_ s - 1 )2^-^+*- s ii Zl >m t -t + s. 

In general, for ^ u C {1, . . . , s} we have 

\{k G X>* : kj > for j G w, = otherwise , fi(k) — Zi}\ 

{ (^itj-'r 1 ) if^<m 4 -t+| M |, 
I (*i+l«l-i)2*i-^+*-l«l if Zl > mi - t + \u\. 

Thus, in general, for z\ > m ; — t + 1 we have 

|{fc G : p(k) = Zl }\ « s ~ ^-m^-i 

Step II Let now k G T>* and z 2 = z — Z\ be fixed and /j, ( ^ ) = i^- for 
1 < j < s - We have j^- — tUj| < 2 and / G X>*/ implies that 

Cl,l/ll,o + • • • + Ci ) 6 1]2 _i/li j 6 li2 _2 + Ci if , 12 + C*1,6 M + 
C*2,l/i2,0 + ' ' • + ^2^2 2-1/12,62,2-2 + £2,62,2 + C2,6 2 ,i + 

C s ,lh s fl + h 0*5^ 2-1^,65,2-2 + C S) 6 s 2 + C S ,6 S]1 = 0, 

where c^o = 0. Note that we consider now to be fixed, thus the hj/s 
are also fixed. Thus it follows that 

Y & ^ + Y + Y &,b j:2 + °i> b i,i) = & > 

for some vector c which is determined by the hij which are fixed since 
k is given. Let fJ,(lj) = Wj for 1 < j < s. Then we have 61,1 + • • - + 6 s ,i = 
Wi H h w s = z 2 . 

Since /ij is fixed by kj for 1 < j < s, it follows that for each given set 
of bj t i, bj j2 , 1 < j < s, at most one such solution exists. 
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Now I G V*, implies that 

61,1 + 61,2 + 62,1 + ^2,2 H h 6 a ,i + fo s,2 > 2m»/ - *, 

^2 = 61,1 + 62,1 H h 6 a ,i > "V - f , 

thus we have 

&i,2 H 1- ^,2 > 2mj/ - i - z 2 + 1. 

Let = 5j + 6^2, thus (5-,- > (where Sj = if lj = 0). Then we have 

Z2 = &i,i H h 6 s ,i 

= (Ji + • • • + 5 S + 61,2 + • • • + & s ,2 

> 81 H h<5 i + 27n i '-*-2 2 + l- 

and therefore 

5i + • • • + 5 S < 2z 2 - 2m v + t. 

Thus the number of bi t i, . . . , b Sj i with + • • • + b Sj i = z 2 is bounded 
from above by 

2z2—2m i /+t . N . 

y~v / /• t ••• - L \ / 1 ! : 2 - 2;/;,' -+ / * 

r=0 



s — 1 / V s 



Now consider the number of possible choices of bj )2 . If j G Ur=-i a ^' 
then 6j 2 is fixed since kj is fixed and if j G a_ 2 , then bj t i > bj 2 > 

Note that 6^3 is fixed since kj is fixed for all 1 < j < s . By the order 
3 digital net property and I G V*, we have 

&1.1 + h,2 + 6l,3 H 1" + & s,2 + & S ,3 > 3771? - t, 

61,1 + 61,2 H 1- + b s ,2 > 2m v - t, 

Z2 = &i,i H h 6 s ,i > mi- - t 

Let 4 = &i j2 + &2,2 H h & s ,2 < -22- Then 

&i,3 H 1- &s,3 > 3mj- - t - 2 2 - 4 > 3mj' - t - 2z 2 - 

Let 6j j2 = 5j + bj j3 , then > 0. Then we have 

Z 2 > 6l,2 H 1" &s,2 
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= 5[ + ■ ■ ■ + 5' s + 6i, 3 + • • • + b S:3 
> 5[ + ■ ■ ■ + 5' s + 3m v -t-2z 2 + l 

and therefore 

5[ + --- + 6' s <3z 2 - 3m;/ + 1 - 1. 
Thus the number of admissible &i i2 , . . . , b s>2 is bounded from above by 

3Z2—3m i l +t 

E 



r + s — 1 \ ( 3z 2 — 3mi> + t + s 
s-1 



r=0 

Thus we have 

(2z 2 - 2m v + t + s\ (3z 2 - 3m v + t + s 

\R(k,z 2 )\ < 

Inserting the results from Step I and Step II into (JSJ) yields 

z-m^+t-l . . 

\JiM\ < £ 2--- + *- i r l +_ s i M 

'2z — 2^i — 2mj/ + t + s\ f 3z — 3z\ — 3m.j/ + t + s 

Let now z = m,j + mj' — 2t + 2 + k for k > 0. Then we have 
1^0^ + "V -2t + 2 + K )| 

C) z> 1 ( z 'i+ m i- t + s 



s-1 



s E 2 

4=0 

2(k-z[) + s + 1- t\ (3{k -z[) + s + 3-2t 
s 



X 



Set 

„, , M n -# fz[+m i -t + s\ /24 + s + l-A /34 + S + 3-2* 

Bfe '^ ) = 2, ( .-1 )( . )( . 

Notice that one can obtain a similar bound with the roles of k and I inter- 
changed. Thus it follows that 

K 

\J i<v (mi + m v -2£ + 2 + k)| < £ min{£(zl, k - z[), B(k - z[, z' x )} 

z[=0 
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\k/2\ 

4=0 

r«/ai 

z[=0 

r*/2i 



< 2K /2+2 /*! + m* - t + s\ 

2(« - 4) + s + 1 - A /3(k - z£) + s + 3 - 2t 



4=0 

x 

We have 6 j 3 > 0^5 for 1 < j < s by Lemma [T] and by the order 5 
digital net property we have (where t(a') is the t-value of the order a' digital 
(t, s)-sequence for 1 < a' < 5) 

a>i,i + 0.1,2 + ai,3 + oi, 4 + oi i5 H h a S) i + a s , 2 + a Si3 + a s , 4 + a s , 5 > 5m; - i(5) 

^1 = a>i,\ H h o Sj i > ai )2 H h o Sj2 > o lj3 H h o Si3 > a M H h a Sj4 . 

Thus we obtain 

z-Zi = z 2 = 6i,iH 1-65,1 > 6i, 3 H h6 s , 3 > oi j5 H ho Sj5 > 5m*-£(5)-4zi 

and therefore 

z > 5m; - t(5) - 3zi > 5m; - 1 (5) - 3(r^ - t(l) + 1) = 2m< - 1(5) + 3t(l) - 3. 
Analogously we have 

z > 2m,/ - £(5) + 3t(l) - 3. 

Thus we have \ Jiy(z)\ = if z < 2max{m;,mj/} — t(5) + 3t(l) — 3. 

Assume that i < z', then since t, t(5),t(l) depends only on the dimension 
s but not on m;, we obtain 

£ kmk. t ^ 

(k,l)ej i>it z=2mj-*(5)+3t(l)-2 

1 I^K + m,, -2t + 2 + K)| 

s Onii+m,/ 



«=m i -m < /+2t-t(S)+3t(l)-4 
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oo r(K+m i -m i ,+2t-t(5)+3t(l)-4)/2] / N 

1 \ 1 \ / Zy + TT^j — t + S 

<s 2 (3m i +m i ,)/2 2^ ^ V S - 1 

^ /2(« - ^) + s + 1 - ^ /"3(k - ^) + s + 3 - 2t' 



Now the sum over z[ can be estimated by the number of terms and by 
substituting the value of z[ in each binomial coefficient which yields the 
largest factor. Thus we obtain 



E KM) I 9(3m . + ™ /)/2 Yl 



J k + rrii — my + 5t 



2(3m i +m i ,)/2 Z-^i 2 K / 2+1 



X 



X 



k + 2rrii — rriii + At + s 
s-1 

2k + s + 1 - t\ [?>k + s + 3 - 2t 



Therefore we obtain 

£lA*s) «s ± E im-mM m ^ )/2 E^ 

l<i<j'<r k=0 

(« + 2mi + 4t + s) 3 " 1 (2k + s + 1 - t) s (3k + s + 3 - 2tf 



x- 



(s-1)! s! 



Now, for fixed i, we have Y?i>=i( m i ~ m v )2^- m ^l 2 < Y^=i<l 2 ~ q/2 < °°- 
Furthermore, we have mi < (log AT)/ (log 2). Therefore 

j r oo 2s+l 

1=1 K=0 
o_] OO o„ 

(log^V)'- 1 

where r = S'(iV) denotes the number of nonzero digits in the binary expansion 
of N. Thus the result follows by taking the square root. □ 
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4 The proof of Corollary U 

To prove Corollary [TJ we first prove a bound on the £2 discrepancy of order 
3 digital nets. 



4.1 A bound on the C% discrepancy of order 3 digital 
nets 

Theorem 2 Let s,m G N. For every (digitally shifted) order 3 digital 
(t,m, s)-netV2 m , s over¥ 2 we have 



£2,2™ (^2™, s) "C* 



??? 



( S -l)/2 



ym—t 



Proof. Let Ci, . . . , C s G F 2 nxm be the generating matrices of P2 m s and recall 
the definition 

V = {keW : Cjh + ■■■ + Cjk s = (mod 2)} 

and £>* = V\ {0}. Set 

J={(fc,/)GD*xF:r(jfc,/)^0} 

and 

J(^) = {(M) e J : fj.(k) + fi(l) = z}. 

Hence J = and J{z) = Ji,i>(z) from the proof of Theorem [I] with i = i' 
and f»i = m. 
We have 



^2 2 m (J^2 m ,s 



r(fc, Z)waL.(cr)walz( 



For (fc, Z) G J" we have 



< E i r ( fc ' z ) 



k,iex>* 



E i r ^) 



|r(fc,Z)| < 



1 



3s2M(fc)+M(0 



Thus we have 



4 2m (p)« s e ^ 



(fc,l)e«7 
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2M(fc)+M0 ' 



It follows from k G T>* that fi(k) > m — t and from I e T>* it follows also 
that > m — t and hence //(fe) + /x(Z) > 2(m — t + 1). 
Thus we have 

ZV(7>)<, £ (10) 

z=2(m-t+l) 

From the proof of Theorem [1] (with % = %' and rrii = m) we obtain 

z—m+t—l 



z—ra-TU—L / \ 



zi=m— 4+1 

'2z — 2^ — 2m + 1 + f3z — 3z x — 3m + t + s 

s 



x 



Let now z = 2m — 2t + 2 + k for k > 0. Then we have 

\J{2m - 2t + 2 + k)\ < B ( z 'n K ~ z 'i)' 

4=0 

where 

2 , fz' x + m - t + (2z' 2 + s + 1 - ^ f 34 + s + 3 - 2i 
s - 1 



B(z[, z' 2 ) = T 
Then 



r«/2i 

1.7(27/1 -2t + 2 + «)| < 2 nrin{S(2i,/c-4),B(«- 4, 4)} 

4=0 

r«/2i 

4=0 



< 2K /2+2 £ / ^ + m - t + s\ 
4=0 ^ 1 ' 



2(k - z[) + s + 1 - A /3(k - 4) + s + 3 - 2t 



Inserting this result into fllOp we obtain 
2 1 ^|J(2m-2t + 2 + K) 

Vl' 2m . s J ^ s 9 2m-2t+2 



22m-2t+2 / / 2 K 

K=0 
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oo . \k/2\ , , 



1 x ^ 1 x ^ ( z[ + m — t + s 

~m-2t / -> Ok/2 / / 



22m— 2t / j 2 K / 2 ' J \ g 1 

k=0 z[=0 



X 



< 



'2(k - 4) + s + 1 - ^ /3(k - z[) + s + 3 - 2f 
1 ^K/2 + l(K+l + m-t + s) s - 1 



2 2m-2t ^ 2 K / 2 (S - 1)! 

(2k + s + 1 - t) s (3/t + s + 3 - 2tf 
X il s! 

Therefore we obtain 

s— 1 00 3s <f— 1 

m x k m 



L 2)2m (P 2 "vJ <s 9 2m-2i W2" ^ 



22m-2t / / 2 K / 2 2 2m " 2t ' 

K = 

Thus the result follows by taking the square root. □ 



4.2 The proof of Corollary QQ 

The proof of Corollary [T] is intimately connected with Theorem [2l 

Proof. For an integer iV > 2 we choose m e N such that 2 m_1 < A^ < 2 m . Let 
P 2 ™,s be an order 3 digital (t, m, s)-net over F 2 with the property that the first 
component of Vi m ,s is a (®i m i l)-net over F 2 . Note that such nets exist for 
every m and can be obtained in the following way: Take the digital sequence 
introduced in Section [L2l in dimension 3s — 1. Concatenate the component 
n2~ m , so that the new point set is of the form {n2~ m 1 y n>1 , y nj2 , . . . , y n ,3s-i), 
where (y n ,i, ■ ■ ■ , Vn,3s-i) is the nth point of the sequence. Apply the digit in- 
terleaving composition to this point set. Then the point set (n2~ m , y n l , y nt2 ) 
for < n < 2 m is a digital (0, m, 3)-net and therefore [TO, Proposition 1] 
shows that the first component of the order 3 digital net obtained this way 
is a digital (0, m, l)-net. 

Note that the quality parameter t only depends on s and not on m. We 
now proceed as in |3]. According to Theorem [2] we have 

m (s-l)/2 

£ 2 ,2«(7V, S ) <s 2m • (11) 
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As shown above, the first component of V 2 m s is a digital (0, m, l)-net over 
F 2 . Hence the subset 

N' 



Vn,s '■— v 2 m yS n 



s-l 



contains exactly N points. We define the point set 



V N , S := 

Then we have 

(NU N (V N , S )) 2 



N ' 



Xi,x 2 , ■ ■ ■ ,x s ) : (x u x 2 , . . . , x s ) G V 



N,s 



\A([0,y),N,V N!S )) - N\ s ([0,y))\ 2 dy 



[0,1] > 



JO 



A[[0, yi N2- m ) xY[[0, yi ),N,V N , s 



i=2 



N 



2 m —yi---y< 



dj/i • • • dy s 



2«i /-7V/2 m /-l 



A([0,y),N,V N , s ) - 2 m X s ([0,y)) 



dy 



N Jo 



^0 



\A([0,y),2 m ,V 2m , s ) -2 m X s ([0,y))\ 2 dy 



2 m 

— -jj (2 m £ 2 ,2 m (7V\s)) • 
With (dlD we obtain 

2 m 

(NC 2 , N (V N , S )) 2 < s <s (logiV) 5 - 1 . 

Taking the square root and dividing by N we finally obtain 

(logiV) (s - 1)/2 



N 



□ 
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